FTUAM 10/95, April 1995 



HEAVY QUARKONIUM SYSTEMS 
AND 

NONPERTURBATIVE FIELD CORRELATORS 



Yu. A. Simonov* 

ITEP, 

B. Cheremushkinskaya 25, 
117259, Moscow, Russia 

S. Titard** and F. J. Yndurain*** 

Departamento de Fisica Teohca, C-XI, 
Universidad Autonoma de Madrid, Canto Blanco, 
E-28049, Madrid, Spain. 



Abstract. 

Bound states of heavy quarks are considered. Using the path integral formalism we are 
able to rederive, in a gauge invariant way, the Leutwyler-Voloshin short distance analysis 
as well as a long distance linear potential. At all distances we describe the states in 
terms of nonperturbative field correlators, and we include radiative corrections at short 
and intermediate distances. For intermediate distance states (particularly bb with n = 2) 
our results improve, qualitatively and quantitatively, standard analyses, thanks mostly to 
being able to take into account the finiteness of the correlation time. 
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Some time ago, Leutwylert-*^! and Voloshint^l showed that the gluon condensate {ag G"^) 
controls the leading nonperturbative effects for heavy qq states at short distances. Their 
analysis was completed in refs. 3,4, in particular by extending it to spin dependent split- 
tings and by including relativistic and one loop radiative corrections, an essential ingredient 
in the analysis. With these additions it was then shown that a consistent description of 
states n — lhh and, to a lesser extent, hh states with n = 2 and cc states with n = 1 could 
be obtained.* 

As already pointed out in refs. 1,2 the approach fails for large n. The reason is that 
nonperturbative contributions grow like n^, quickly getting out of hand. For example, for 
the spin-independent spectrum we have. 



log + ^(n + / + 1) 



M(n,0 = 2m{l 



+ 



87rn2 2(mCF5s)- 



}. (I) 



Here I is the angular momentum, m is the pole mass of the quark, Po = (33 — 2n/)/3, 
Ci? = I and 5s embodies part of the radiative corrections: 

93 — lOn/ 7e/3o 



Here is the renormalization point. Finally, in the leading nonperturbative approximation, 
and assuming a constant gluon condensate density, 

the Sni are numbers of order unity: 

_ 624 _ 1051 _ 769 456 _ 9 929 

tlO 425' 20 ) t30 463 239' 21 9 944'" 

(2) 

We will not consider in this note relativistic corrections. 

Clearly, the nonperturbative correction in Eq.(l) blows up very quickly for increasing 
n making the method totally unsuitable already for bb with n = 3, and cc with n = 2. 
The way out of this difficulty found in most of the literature is to use a phenomenological 
potential, 

err, (7^ ~ 0.45 MeV. (3) 

Now, and although a linear potential yields a correct description of long distance qq forces, 
the methods lack rigour in that, as is well known t^'-*^'^!, a linear potential is incompatible 
with known QCD results at short distances, where indeed it does not represent a good 
approximation of e.g., the nonperturbative part of (1). 



* n is the principal quantum number. We will use standard atomic spectroscopic 
notation 
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Because of this it is desirable to develop a framework which, in suitable limits, implies 
both the Leutwyler-Voloshin short distance results as well as the long distance description 
in terms of a linear potential. This framework is an elaboration of that developed in refs. 7, 
where it was shown from first principles how one can derive a linear potential from first 
principles. In the present note we explore the short and intermediate distances, where we 
rederive the Leutwyler-Voloshin description, improved both by getting better agreement 
with experiment for the states where it is valid, and extending its range of applicability 
to intermediate distances. The reason for this improvement lies in that our treatment 
includes the nonlocal character of the gluonic condensate in the form of a finite correlation 
time, Tg. The results of refs. 1 to 4 are then recovered in the limit — > oo. 

The nonlocal condensates have been considered previouslyt^'"'^^^, in particular with 
the aim of finding this correlation time. In this note we improve upon the treatment of 
refs. 9,10 first by using a Lorentz and gauge invariant path integral formulation which 
would allow us, if so wished, to incorporate relativistic and spin effects**. Secondly, we 
incorporate radiative corrections which are essential to give a meaning to parameters like 
as{fi) or m. 



Description of the method. 

The method uses the path integral formalism. Because in this note we are only 
interested in nonrclativistic, spin independent splittings, we start directly with the nonrel- 
ativistic qq Green's function'^'^L For large time T, 

G{x,x-y,y) =4:7X1^6-^"^^ J VzVze-^^°+^°^ {W{C)), (4) 

W{C) is the Wilson loop operator corresponding to the closed contour C which includes the 
q, q paths. W{C) should also include initial and final parallel transporters, ^{x, x), y) 
with e.g., 

^{x,x) = P expi^ / dzi^Bi^{z), (5) 

J X 

P denoting path ordering. Actually, we can omit the parallel transporters, and at the same 
time avoid problems with the renormalization of the Wilson loop by choosing x — x, y = y, 
which will prove sufficient for our purposes. 

In order to take into account the nonperturbative character of the interaction we split 
the gluonic field B as 

Bi^ = b^ + a^. 

** The methods to acomplish this would be like the ones developed in the fifth paper of 
ref. 7 and in ref. 8 
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The separation will be such that, by definition, the vacuum expectation value of Wick 
ordered products of vanishes, so the correlator may be written in terms solely of b^: 

{G{x)G{y)) ^ {G,{x)G,{y)), (6) 

and Gb is constructed with only the piece of B^. One may expand in powers of the 6^ 
and thus write the Wilson loop average as 

{W{C)) = J VaPexpig J^dz^a^ 

= Wo + W2 + ..., (7) 

where the transporter $a is given by an equation like (5), but in terms of only the piece 
of B^. 

Let us discuss the first term in the r.h.s. of Eq. (7). Using the cluster expansion we 

get 

Wo = Zex.p{(j)2 + (f)4 + ...), 



and regularization (to be absorbed in Z) is implied in this integral. It should be noted that 
(j)2 contains all ladder- type exchanges, and in addition also "Abelian crossed" diagrams - 
those where the times of the vertices can be not ordered, but where the color generators 
t?^ are always kept in the same order. Because of this, all crossed diagrams (with the 
exception of the "Abelian crossed") are contained in ^4. It is remarkable that each term 
(f)2n in Eq.(8) sums up an infinite series of diagrams. In particular, and as we will see below, 
exp(/)2 contains all powers of as/v {v being the velocity of the quarks) so the calculation 
is exact in the nonrelativistic limit. 

For heavy (and slow) quarks, (8) becomes, 

i.e., a singlet one-gluon exchange potential, as expected. 

We then turn to W2. When expanding it in powers of one gets typical terms like 

C% Tr(t«i . . . t''" e hi t"" . . . 6^ ^^ . .). 

Because of the equality 

tH^'t" = -t''/2Nc, 
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one obtains, for all exchanges in the time interval between the times of 6^(2;) and bi,{z') a 
factor —1/Nc instead of the factor Cp that is found for other exchanges. As a result we 
may write W2 as 

W2 = ^ f dz„, I dz', 



c 



'C 



X exp 



rdtv^'^+ rdtv^'^+ f^dtv^'^ 

Jt2 Jo Jtl 



(10) 



with V^Q^^ = CpOis /r, V^Q^^ = —{1/Nc) cx-s /f the singlet, octet potentials respectively. The 
relevance of the octet potential was already noted in refs. 1,2. It appears in our derivation 
in a fully gauge invariant way, in connection with a gauge invariant Green's function with 
the gauge invariant quantity W2 as the kernel. 



Evaluation of W2 

One first uses the Fock-Schwinger gauge (see ref. 6 for details, including a modification 
of this gauge) to write, 

/ dz^ / dz'^{b^b^)= / dap_pda'^y^{Gi,^p{z)Gi,^x{z')), 
Jc Jc J 

and the da are surface differentials. Including also parallel transporters, equal to unity in 
the F-S gauge, we find the gauge-invariant expression, 

j dz^, J dzl {bi^b^) = J dap_pdal^l^^{xo,w)Gb,^p^w,xo) 

x^xo,w')Gb„x^w',xo)y (11) 

As was shown in ref. (7) xq may be chosen between w and w', up to additional contributions 
of order 6^, that we are neglecting here. Then we divide the total time interval T into 
three parts (c/. Fig.l): 

(/) 0<t<w'4^ 

(II) wi<t<W4 

(III) W4<t<T. 

Separating out the trivial cm. motion we get, in regions (i), {III), the singlet Coulomb 
Green's function. 
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Fig. 1 -The surface C 



G^^)(r(ti),r(t2);ti -^2) = / Vr{t) exp 



-m 



/ dtr'^ + Cpas I 

Jt2 Jt2 



dt 

t2 W) 



for quarks of equal mass (so that the reduced mass is m/2). In region (//), however, 
we find the octet Green's function, 



G[?^(r(w4),r(w4);w4 - W4) = / Dr(t) exp 



Jw' 



rW4, 
Jw' 



dt 



Inserting now this into (10) yields the correction to the total Green's function, G = G'^^^ + 
SG: 



SG = 



J d^r{w4)G^'^\r{T),r{w4);T - W4) J d^r{w'4)G^^\r{w4),r{w'4) ■,W4 - W4) 
da^,{w) J dapx{w'){G^,{w)Gpx{w'))G^^\r{w'4),r{0y,w'4). (12) 



At this point it is convenient to specify the surface inside contour C, which we do by 
connecting the points z{t) and z{t) by a straight line. In the cm. system, 

dcr^i/ = a^^dzodp, 

and, in the nonrelativistic approximation, the a^- may be neglected. Thus, and as expected, 
only the chromoelectric piece £ of G^j^ survives in the correlator in Eq. (12). We then 
expand this correlator in invariants: 

{g^£i{x)Sj{y)) = ^[%A(a; - y) + h^hj dD^/dh% 
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(13) 

h = x-y, 
and 



/^{z)^D{z) + D,{z)+z'^. (14) 
The invariants D, Di are normalized so that 

L'(0) + Di(0) = 27r(a,G2). (15) 
Inserting then (13) into (12) gives, 

6G = ^ J d^r J dVy udp j r^d/3'G^^^(r(T),r) GP{r,r')G^c \'^' ^'^i^))- 

We next consider the matrix element of 5G between Coulombic states, \nl). The resulting 
expression simplifies if we use the spectral decomposition for Gc'- 

G^c'^\ry;t) = {r\eM-H^c'^h]\r') 

k 

Identifying the energy shifts from the relation 

G = G(^) + 5G ~ G^^^(l -T5Eni) 

valid for T ^ oo we get 

9^ / d/5d/?'A(p)5](nZ|r,e-(/^-^)r|/c(8)) 
x ^(3) l^S) . (fe(8)|r-e^P(^-^)'-V0- (16) 

Ef, — En — ip4 

The states |A;(8)) are eigenstates of the octet Hamiltonian, 

= p2/TO + as /2Nc r 

with eigevalues Ej^\ ^{p) is the Fourier transform of A (a;). 

Eq. (16) is our basic equation. The correlator A(a;) depends on x as 

A{x) ^ f{\x\/T,) 

and is expected to decrease exponentially for large |a;|, in Euclidean space. The correlation 
length Tg may be related to the string tension, 

T,-' - = ^ ^^^^ - 0-35 GeV. (17) 
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(for the derivation of this equation, see below). We have now two regions. For very heavy 
quarks, and small n, 

Then we approximate A(a;) ~ constant, and hence A(p) ~ ^4(p) so that 

_ TTjasG^) ^ {nl\ri\k{8)) {k{8)\n\nl) 

^^nl - 2^ (8) (5) 

This is the equation obtained in refs. 1,2 and, upon calculating the r.h.s. of (18) one 
indeed finds the nonperturbative piece of Eq. (1). 

The improvement over (18) represented by (16) lies in that it involves the correlator 
{Ei{x)Sj{y)), i.e., one takes account of nonlocality of the correlator. To implement this it 
is convenient to distinguish two regimes: i) We consider that fiT is smaller than 1/a = 
mCpOis/'^i but one has /xt — \En\- ii) jJix 3> \En\. 

We will first consider case (i). Then one can neglect |p| in the exponents of Eq. (16) 
but ht should be kept in the denominator there. Using now an exponential form for A (a;) 
we obtain. 

Substituting this into (16) we get the energy shifts, 

7r(a^G^ rAnl). (20) 

This is precisely the approximation postulated in refs. 9,10. Clearly, as Tg — > oo {ht — > 0), 
(20) reproduces (18); but, as we will see, (20) represents an important improvement over 
(18) both from the conceptual and the phenomenological point of view in the intermediate 
distance region. 

Then we turn to the regime {ii), (it ^ \Eif^\. In this case the velocity tends to zero, 
the nonlocality of the interaction tends to zero as compared to the quark rotation period 
(which in the Coulombic approximation would be Tg = 1 / 1 En | ) , and the interaction may 

therefore be described by a local potential. In fact: considering Eq. (16), it now turns out 

( s) . . . (&) 

that we can neglect both En and the kinetic energy term in E^, (indeed, all of it) as 

compared to —ip4. Then one gets 

-j^ /•ri{w4) /•ri(w'^) POO 

5Eni ^ —{nl\ / dwi / dw[ / d{w4 - w'4^)A{w - w')\nl) 
36 Jo Jo Jo 

and in the limit in which we are now working we can approximate ri{w'^) ~ rj(w4). One 
obtains here the matrix elements of a local potential which may be written in terms of 
D, i?i:m 



t^(r)= 3^{2r dA ^ duD{X,u) 
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pr 

+ 



pr poo 

J AdA y diy[-2D{X,u) + Di{X,u)]j. 



Note that both Z), Di depend on A, v through the combination X'^ + u'^. At large r, and 
as this equation shows, U (r) behaves like 

U{r) ~ err — Constant, ^ ~ ^2 J dicida;2 -D(a;i, a;2). 

If we now used the ansatz (19), we would obtain the announced relation between ht and 
a. Of course in this situation the strategy of treating the effects of the gluon condensate 
as a perturbation of the Coulombic potential is no more appropriate. One should rather 
take U{r), together with the Coulombic potential, as part of the unperturbed Schrodinger 
equation. We leave the subject here referring to the various existing analyses t^'^-*^'^^! for 
details. (In particular, refs. 11 are the ones closer in spirit to the work here in what 
regards the treatment of the nonperturbative effects, while ref. 12 incorporates radiative 
corrections to a phenomenological long distance potential) 



Phenomenology 

For the phenomenological analysis we will generalize Eq. (1) by writing. 



M{n,l) = 2m|l 



8n2 



2 



log +i/;{n + l + l) 



87rn2 2{mCFOisY 



where the eniipr) are obtained solving Eq. (20); for ht —>■ 0, eni{0) = €ni, this last being 
the quantities given in Eqs. (1), (2). 

Unlike in the case //t = 0, where a closed expression could be found for the eni, the 
^nlilJ'T) may only be computed numerically. However, a fairly precise evaluation may be 
obtained by neglecting the potential /2Ncr in (20), then working with p-space Coulomb 
functions. To a very tolerable ~5% accuracy it follows that we may approximate, 

eni _ / 2n \ ^ iJ,T 



enlM - , r/n ^ — ^ ^, (22) 

where, 

pio = 0.62, P20 = 0.76, P30 « 0.9, p2i = 0.70. 
Substituting into Eq. (21) then gives us a very explicit generalization of Eq. (1), valid in 

(S)i 

the intermediate region /it ~ \En \- 



9 



For the numerical calculation we proceed as follows. For bb we take the optimum 
values for the renormalization point given in ref. 4 (which were obtained neglecting //t)- 
We therefore choose, 



For mixed n we take the value corresponding to the smaller n. For cc we, somewhat 
arbitrarily, choose /j, = 0.95 GeV. For the basic QCD parameters we take 



and thus the corresponding values of are as{l.5 GeV) = 0.27, o;s(0.95 GeV) — 0.35. 

We will not consider varying these quantities. A variation of A can be largely com- 
pensated by a corresponding variation of /x and likewise, and because the e„i(//T) depend 
almost exactly on the ratio {as C^) / ht-, a variation of the condensate may be balanced by 
a compensating variation of the correlation time, Tg = /j,^^. 

We now have two possibilities: fit fix to each individual splitting, and compare the 
resuls among themselves and with the one coming from the string tension; or take fiT from 
the last, Eq. (17) and then predict the splittings. If we do the first we find 

liT = 0.40 GeV (2S - IS), /it = 0.76 GeV (3S - 2S), /it = 0.59 GeV (2S - 2P). (23) 

For the cc case we only consider the 2S-1S splitting, where we get 



Clearly, the more reliable calculation is that of the 2S-1S bb splitting: not only the radiative 
corrections are known (unlike for the 2S-2P case) but also it falls inside the conditions of 
regime {i), unlike the 3S-2S splitting and, even more, the 2S-1S cc one. It is then gratifying 
that the value of /it that follows from the bb 2S-1S splitting, /jt = 0.4 GeV, is the one 
which is in better agreement with the value /it = 0.32 GeV obtained with the comparison 
with the linear potential, Eq. (17). 

If we now choose the second possibility, we fix fj,T- The corresponding results are 
summarized in Table I. 



11= 1.5GeV,forn= 1; // = 0.95 GeV, forn = 2. 



A(n/ = 4, 2 loops) = 200 MeV, {asG^) = 0.042 GeV, 



liT = 1.23 GeV. 



(24) 



splitting 



I^T = iiT = 0.32 //T = 0.40 



exp. 



2S - IS (bb) 
2S - 2P (bb) 
3S - 2S (bb) 



479" 
181" 
4 570 



590 
162 
748 



522 
147 
614 



558 MeV 
123 MeV 
332 MeV 



2S - IS (cc) 



9 733 



1930 



1626 



670 MeV 



Table I.- Predicted splittings, and experiment, 
(a): Values from ref 4, with /j, ~ 0.95 for both 
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The n = 2 bb states are certainly better described than with the approximation Tq = 
oo of refs. 1, 2, 4. Particularly important is the fact that inclusion of the finite correlation 
time stabilizes the calculation. For example, if we had taken ^ = 1.5, and fix = 0, for 
the 2S-1S splitting for bottomium, we would have obtained the absurd value of 1 944 
Ge V. In what respects the n = 3 and the n = 2cc states the improvement is marginal, 
in the sense that the basic assumption, vu., that one can treat the nonperturbative effects 
at leading order fails, as is obvious from the figures in the column " fiT = 0" in Table I. 
Indeed these states fall clearly in regime (ii) and should therefore be better described with 
a local potential as discussed extensively in the existing literature, of which we, and for 
illustrative purposes, single out ref. 11, where the nonperturbative effects (including spin- 
dependent splittings) are treated with methods like ours, but where radiative corrections 
are ignored; or ref. 12 where radiative corrections are incorporated but the confining 
potential is introduced phenomenologically. 

We would like to end this note with a few words on extensions of this work. An 
obvious one is to include the treatment of spin effects, and a calculation of the wave 
functions. Then it would be very desirable to evaluate the radiative corrections to the 
nonperturbative terms, as this would greatly diminish the dependence of these terms on the 
renormalization point, /U, thereby substantially increasing the stability of the calculation. 
Finally, and to be able to extend the calculation to intermediate distances with success, 
one should abandon the treatment of nonperturbative effects at first order: an iterative 
approach should certainly yield better results. 
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